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ABSTRACT 

Low-mass starbursting galaxies have been proposed as the explanation of the excess 
of faint galaxies observed at intermediate redshifts. If this hypothesis is correct, then 
intermediate redshift galaxies should rotate more slowly than nearby galaxies with the 
same rest-frame luminosity. We present the results of a survey of the internal kinematics 
of intermediate redshift (z = 0.25—0.45) field galaxies to search for this effect. Using the 
Canada-France-Hawaii Telescope, spatially-resolved spectra of the [O II] AA 3726—3729 
A doublet emission line have been obtained for 22 galaxies. V ro t sin i and [O II] disk 
scale lengths have been extracted from each galaxy spectrum using a Bayesian fitting 
technique. 

Galaxies in the sample are found to be ~1.5— 2.0 mag brighter than expected from 
their rotation velocity and the local Tully-Fisher (TF) relation. Low-mass galaxies ex- 
hibit a wider range of evolution relative to the TF relation than high-mass galaxies. The 
main source of uncertainty in this result is the large scatter in the local TF relation for 
late-type galaxies. Luminosity-dependent luminosity evolution neatly reconciles the lack 
of evolution seen in other works with the results of our survey. It is also found that the 
overall properties of [Oil] kinematics at intermediate redshifts are varied. For example, 
25% of the field galaxies in the sample have [Oil] kinematics unrelated to rotation; [Oil] 
emission is confined to the nucleus in most of these galaxies. Anomalous kinematics 
is found to be related to the presence of companions — i.e. minor merger events. A 
Doppler ellipse similar to those found in local dwarf galaxies has been observed in a z = 
0.35 galaxy, and may be interpreted as a supernova-driven supershell. 

Subject headings: galaxies: kinematics and dynamics — galaxies: evolution — methods: 
data analysis 
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1. Introduction 

The observed number counts of faint galaxies ex- 
ceed the predictions of standard models in which evo- 
lutionary effects are ignored. Yet the faint galaxy red- 
shift distribution appears to be well modelled by these 
same no-evolution models (Lilly et al. 1995, Lilly 
1993, Koo and Kron 1992, and references therein). 
The amount of evolution that faint galaxies have un- 
dergone (relative to galaxies at the present epoch) 
is clearly a critical parameter for understanding the 
Universe at intermediate redshifts. Here we focus on 
the novel technique of internal kinematics for study- 
ing galaxy evolution. 

Internal kinematics is directly related to a sin- 
gle fundamental property of galaxies: mass. Obser- 
vations of the velocity field of intermediate rcdshift 
galaxies can therefore be used to quantify the amount 
of luminosity evolution that these galaxies have un- 
dergone. As an illustration, consider a galaxy with 
rest-frame absolute luminosity equal to that of the 
Milky Way. If this galaxy were also as massive as the 
Milky Way (i.e. if it were unevolved in luminosity rel- 
ative to the Milky Way), then a rotation velocity of 
~ 200 km/s would be observed. On the other hand, 
if it were in fact a lower mass object that had been 
boosted in luminosity by ~10 x (as suggested in the 
luminosity-dependent luminosity evolution scenario - 
Broadhurst et al. 1988), then the TF relation would 
predict a rotation velocity ~ 100 km/s. 

This kincmatical approach to measuring evolution 
is direct. It is not affected by uncertainties in mod- 
els based on local luminosity functions. Furthermore, 
whereas luminosity functions derived from rcdshift 
surveys constrain the average amount of evolution in 
an entire galaxy population, the present approach can 
measure luminosity enhancement in individual galax- 
ies, and can thus tie luminosity evolution directly to 
galaxy properties on a galaxy by galaxy basis. 

Most previous kinematical studies (Franx 1993, 
Vogt et al. 1993, Colless 1994, Koo et al. 1995, 
Forbes et al. 1996) have suffered either from small 
sample size, or from a lack of spatial information; 
but they have clearly demonstrated that current tele- 
scopes are adequate to the challenge of measuring the 
internal kinematics of intermediate redshift galaxies. 
Recently, Vogt et al. (1996) presented a beautiful 
set of rotation curves (observed with the Keck Tele- 
scope) for nine faint field galaxies in the redshift range 
0.1 < z < 1. These rotation curves appeared similar 



to those of local galaxies in both form and ampli- 
tude. The kinematics of the Keck galaxies showed ev- 
idence for only a modest increase in luminosity (AM B 
< 0.6 mag) compared to the local Tully-Fisher rela- 
tion. This is in apparent disagreement with the strong 
evolution in galactic disks observed by Schade et al. 
(1996) over the range 0.1 < z < 0.6. 

In this paper, we present the internal kinemat- 
ics of a sample of intermediate redshift galaxies, 
and demonstrate how luminosity-dependent luminos- 
ity evolution can reconcile observations of surface 
brightness and internal kinematics. 

2. Observations 

Twenty field galaxies and two cluster galaxies in 
the redshift range 0.25—0.45 were selected from the 
CNOC cluster survey database (Yee et al. 1996). 
This database contains position, Gunn r magnitude, 
g — r color, and redshift for thousands of galaxies. 
Elongated objects were preferentially selected to min- 
imize sin i effects. No size/color selection criteria were 
used. Rest-frame [Oil] equivalent widths W3727 were 
chosen to lie between 20 and 50 A. Target galaxies 
may be [Oil] strong relative to local early-type (Sb 
and earlier) spirals, but they have the same W3727 as 
many local late-type (Sc and later) galaxies (Kenni- 
cutt 1992). Our range of W3727 is also representa- 
tive of the excess galaxy population at intermediate 
redshifts (Broadhurst et al. 1988, Broadhurst et al. 
1992). 

Spatially-resolved observations of the [O II] AA 
3726.1, 3728.8 A doublet emission line were obtained 
with the Multi-Object Spectrograph (MOS) and the 
Subarcsecond Imaging Spectrograph (SIS) at the Ca- 
nada-France-Hawaii 3.6-m Telescope (CFHT) in July- 
August 1994. SIS was ideally suited for studying [Oil] 
kinematics because SIS "tip-tilt" corrections yielded a 
seeing FWHM of 0".5-0".6. We used the B600 grism 
and a LORAL 2048x2048 CCD (QE = 23%, N fl = 8 
e~ ) . The detector was binned by 2 along both axes to 
improve S/N ratio. Image scale and dispersion were 
0".31/pixel and 1.58 A/pixel (MOS) and 0".17/pixel 
and 0.88 A/pixel (SIS). Typical total integration time 
per galaxy was 7200 seconds split in 2—3 exposures 
for cosmic ray removal. 

3. Synthetic Rotation Curve Fitting 

Flux levels in [O II] were very low (typical S/N 
ratio per pixel ~ 2—3); it was therefore important 
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to choose a method of analysis which used all of the 
pixels simultaneously to statistically find the best pa- 
rameter values and their respective uncertainties. We 
adopted a synthetic rotation curve fitting technique, 
in which models were compared to the two dimen- 
sional distribution of pixel intensities in the [O II] line 
(the two dimensions being wavelength or velocity, and 
position along the slit) . The parameters of the fitting 
model were the projected rotation velocity V rot sin 
i of the galaxian disk in km/s, the [Oil] exponential 
disk scale length id in h _1 kpc, the [Oil] total line 
flux in DU, and the dimensionless [Oil] doublet in- 
tensity ratio 1(3726 A)/I(3729 A). (It was necessary 
to include this doublet intensity ratio in the model 
because the [O II] doublet is partially resolved due to 
(1+z) spectral stretching.) 

The 2D model distributions of pixel intensities 
were constructed assuming that the [Oil] emission 
was distributed in a thin exponential disk with a flat 
rotation curve. These [Oil] model disks were con- 
volved with a point-spread-function (PSF) extracted 
from direct images; the effect of placing a slit in front 
of the image was then computed. Finally, the flux 
passing through the synthetic slit was convolved with 
the spectrograph instrumental line shape, which was 
extracted from comparison arc lines. 

The best fitting parameter values were found us- 
ing the Metropolis algorithm (Saha and Williams 
1994). This Bayesian algorithm starts with some set 
of parameter values u> and the associated P(co\D, M) 

the probability that u is the true parameter set 
given the data D and a model M. It then picks a 
possible change 5u> in the parameters and computes 
P{uj + 5lu\D, M). If P{uj + 5u\D, M) > P{us\D, M), 
then the change 5ui is accepted. If P{u> + Slu\D, M) < 
P(oj\D, M), then the change Sco is accepted only some 
fraction P(u)+5lj\D, M) / P(lj\D, M) of the time. Af- 
ter hundreds of iterations, the distribution of accepted 
ijj values will converge to P(w\D, M) provided that all 
possible u> are eventually accessible. Parameter space 
is thus sampled with a density proportional to the 
likchood. 

The trial changes 5u are chosen at random. If they 
are too small (i.e. the parameter search is too "cold"), 
then all iterations will accept changes. If the search 
is too "hot", then none of the iterations will accept 
changes. The "temperature" of the search is regu- 
lated so that half the iterations accept changes. As 
the algorithm proceeds, trial values are stored to si- 
multaneously derive the best parameter values and 



their respective error distribution. For each param- 
eter, we took the median value of P(co\D,M) as the 
best value and the 16th and 84th percentile values 
as our "1 er" error bars. We tested confidence inter- 
vals and parameter value recovery using a set of 342 
simulated rotation curves. We chose input parameter 
values in these simulations to yield worse S/N ratios 
than seen in the data. No biases were detected. The 
Metropolis parameter estimates were also more ro- 
bust to noise than least-squares estimates at low S /N 
ratios. 

4. Results 

4.1. [Oil] Morphologies 

[Oil] kinematics at intermediate redshifts is varied. 
Seven galaxies have [O II] scale lengths significantly 
smaller than their broad-band scale lengths, consis- 
tent with central unresolved [O II] sources. Two of 
these "unresolved [O II]" galaxies are serendipitously 
observed cluster galaxies (we cut extra slitlets on the 
mask for cluster galaxies whenever possible). Such 
discrepancies indicate that the [Oil] gas kinematics 
is decoupled from galaxy rotation in some galaxies 
and could be confined to the nucleus. Two properties 
of kinematically anomalous galaxies have emerged so 
far: (1) All but one have close companions. This sug- 
gests that enhanced star formation activity may be 
the result of merger events. (2) Some of these galax- 
ies appear to be of early- type. Kinematically "anoma- 
lous" field galaxies make up 25% of the field sample. 
This is similar to the fraction of blue-nucleated galax- 
ies observed in HST images of galaxies at z ~ 0.6 
(Schade et al. 1995). 

In one galaxy at z = 0.35, the [Oil] line was donut- 
shaped. This line shape could be interpreted in two 
ways: either (1) the line is made of two rotation 
curves from two galaxies very close together, or (2) 
the line arises from an expanding supershcll, presum- 
ably driven by supernova winds from a massive star- 
burst. If option (2) is correct, then the supershcll has 
a diameter of 2.5 h _1 kpc and a rest- frame expansion 
velocity of about 80 km/s. These characteristics are 
amazingly similar to those of supcrshells observed in 
local dwarf irregular galaxies (Marlowe et al. 1995, 
Martin 1995). In another galaxy at z = 0.42, there is 
a relatively strong [Oil] source unresolved both spec- 
trally and spatially superimposed over a fainter ro- 
tational [O II] component. The source could be a 
giant HII region 2.6 h _1 kpc from the center of the 
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galaxy. More details on these unusual sources will 
appear elsewhere (Simard and Pritchet 1996). 

4.2. An Intermediate Redshift Tully-Fisher 
Relation 

Figure [l] shows kinematical evidence for luminosity 
evolution at intermediate rcdshifts. The locus of the 
local H Q rotation velocity— B band Tully-Fisher (TF) 
relation for all morphological types (from the data of 
Mathewson et al. 1992) is used as a reference. Solid 
circles are the rest-frame V rot sin i 's versus rest-frame 
B magnitude for the kinematically normal galaxies in 
our sample — i.e. galaxies with [Oil] scale length 
consistent with their broad- band scale length. B- 
band k-corrections were computed using the tables of 
Frei and Gunn (1994). The upper long dashed line 
is an unweighted linear fit to all local morphological 
types. This linear fit was then shifted by AMb = 
— 1.0 mag (middle dashed line) and AMg = —2.0 
mag (lower dashed line) to represent various degrees 
of luminosity evolution. 

Galaxies in our sample are ~1.5— 2.0 magnitudes 
brighter than expected from their rotation velocity, 
given the local TF relation defined by the Mathew- 
son et al. sample. Such a large local sample clearly 
shows the scatter in the local TF relation. Locally, 
there are no systematic shifts between loci of differ- 
ent morphological types, but it is obvious that certain 
types (e.g. T = 6) display larger dispersions. 

5. Discussion 

The scatter in the local TF relation has a direct 
impact on the exact amount of luminosity evolution 
derived from our observed TF relation at intermedi- 
ate z. This is because scatter in the local TF re- 
lation may correlate with star formation rates, and 
hence emission line strengths. In order to accurately 
measure magnitude offsets relative to the local Tully- 
Fisher relation, each galaxy in our CFHT sample (or, 
for that matter, any other sample) should be com- 
pared to the local Tully-Fisher relation for galaxies 
with similar [Oil] emission line strengths. This is not 
possible at present. A large, homogeneous sample of 
local emission-line strengths, rotation velocities, mor- 
phologies and absolute magnitudes is required to set- 
tle this issue. This will remain a major limitation of 
internal kinematics studies at high rcdshifts as long as 
technical requirements restrict such studies to strong 
emission-line objects. 



The luminosity-dependent luminosity evolution sce- 
nario neatly reconciles the various amounts of lumi- 
nosity evolution seen in the surface brightness and 
internal kinematics studies. At the low end of the 
galaxy mass spectrum, Koo et al (1995) looked at 
compact, unresolved galaxies with narrow emission 
line widths. The line widths of these galaxies were 
2— 3x lower than expected from the TF relation of 
normal spiral galaxies, implying large amounts of 
evolution (~ 3 mag). These low-mass galaxies had 
lincwidths similar to those of local HII galaxies. At 
the high end of the galaxy mass spectrum, the Keck 
study of Vogt et al. (1996) selected galaxies that were 
intrinsically large (r^ > 3.0 kpc) and bright (Mb < 
—20.7). These galaxies were also intrinsically massive, 
with typical rotation velocities of 200 km/s. The in- 
crease in B luminosity with respect to the local TF 
relation was less than 0.6 mag. 

Our CFHT sample occupies a niche in size and 
mass between the above two Keck samples. The 
CFHT galaxies are typically a full magnitude fainter 
than the objects in the Vogt et al. study. They 
are also intrinsically smaller with typical disk scale 
lengths less than 2.0 kpc, and some of them are barely 
resolved. They are also less massive, with rotation 
velocities ~ 100 km/s, but they are nonetheless more 
massive than the galaxies observed by Koo et al. If 
mass is taken as an indicator of the luminosity that a 
galaxy would have had in a quiescent phase, then all 
three internal kinematics studies can be understood 
with luminosity-dependent luminosity evolution. 

The B-band surface brightness /io(B) of field disk 
galaxies undergoes strong evolution over the redshift 
range 0.1 < z < 0.6 compared to two different local 
log r,j— Mab(B) relations: the Freeman law and an 
empirical z = 0.06 relation for galaxies in Abell 2256 
(Schade et al. 1996). At a redshift of 0.5, A/x (B) is 
equal to —1.1 mag. This is consistent with or perhaps 
slightly less than the evolution seen in our CFHT sam- 
ple, but it is certainly more than the amount of evo- 
lution seen in the Keck sample of Vogt et al.. Looking 
at Figure 1 of Schade et at, there is a hint that sur- 
face brightness evolution depends on disk scale length: 
smaller galaxies evolve more drastically than large 
galaxies. The effect is particularly noticeable in the 
highest redshift bin where the log r^— M^b(B) rela- 
tion clearly curves "down". Taking Mg ~ —21 and 
<Td> = 4.3 kpc (H = 75), one can see on the Schade 
diagram for 0.45 < z < 0.65 that a number of bright, 
large galaxies at (log id = 0.8, M^^B) = —21) show 
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little or no evolution - similar to what is observed in 
the Keck sample! 

The diversity of [Oil] kinematics and the signifi- 
cant amount of luminosity evolution seen in our sur- 
vey illustrate the power of internal kinematics as a 
probe of the nature of intermediate redshift galax- 
ies. Although samples are small compared to those of 
redshift surveys, they already provide exciting direct 
evidence that luminosity-dependent luminosity evo- 
lution may indeed be at the root of the faint galaxy 
excess problem. Tying internal kinematics to other 
galaxy properties on a galaxy by galaxy basis will play 
a key role in our understanding of galaxy evolution at 
intermediate redshifts. 

We thank the Canadian Time Allocation Com- 
mittee for the CFHT for generous grants of observ- 
ing time, and the CFHT organization for the tech- 
nical support which made this project possible. We 
also thank the Canadian Astronomical Data Center 
(CADC) for generous allocations of computer time. 
We gratefully acknowledge financial support from the 
Natural Sciences and Engineering Research Council 
of Canada through an operating grant to CJP and a 
Postgraduate Fellowship to LS. 
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Table 1 

Target galaxies: Characteristics and Results. 



ID z Gunn Gunn W3727 Mg V ro t sin i r^ 

r (5 - r) obs (A) (km/s) (h" 1 kpc) 

(1) (2) (3) (4) (5) (6) (7) (8) 



A2390-101033 ac 


0.2460 


18.90 


0.70 


55 


-20.3 


-212 + 9 

— 11 


0.09 


A2390-100686 11 


0.3822 


21.22 


0.67 


21 


-19.4 


91 + f Q 

— 39 


2.3 


A2390-350416 11 


0.2558 


20.14 


0.27 


35 


-19.5 


62 t 11 


1.9 


A2390-350471 11 


0.2559 


19.40 


0.28 


21 


-20.3 


171 t it 


2.2 


E1512-301037A 11 


0.3457 


21.41 


0.19 


41 


-19.1 


77 + 22 
' ' - 24 


1.8 


E1512-301037B n 


0.3457 


21.41 


0.19 


41 


-19.1 


86 t J 


0.8 


E1512-101526 a 


0.4026 


20.36 


1.12 


48 


-20.2 


999 +143 
zzz - 33 


0.1 


E1512-201429 n 


0.4231 


20.44 


0.56 


26 


-20.4 


151 ± » 


1.0 


E1621-100515 


0.3455 


20.00 


0.42 


22 


-20.3 


-237 + 18 7 


0.8 


A2390-100225 


0.3829 


21.65 


0.89 


26 


-18.9 


61 + 121 

01 -107 


0.6 


A2390-101084 ac 


0.2302 


17.31 


0.68 


110 


-21.7 


-271 t \l 


0.4 


A2390-200928 n 


0.2645 


21.51 


0.22 


40 


-18.2 


65 t\l 


0.3 


A2390-200802 11 


0.3208 


21.16 


0.37 


25 


-19.1 


-198 ± It 


1.0 


A2390-200372 n 


0.3485 


20.15 


0.65 


16 


-20.2 


-86 ± It 


3.2 


E1512-201845 a 


0.3387 


19.60 


0.82 


36 


-20.5 


385 + 11 


0.4 


E1512-201773 


0.3383 


20.06 


0.35 


26 


-20.3 


-116 + {I 


0.9 


E1512-200730 a 


0.4266 


21.35 


0.58 


32 


-19.5 


-H 5 


0.2 


E1512-200334 n 


0.4142 


21.81 


0.70 


36 


-19.0 


60 1 SS 


1.5 


E1512-200672 a 


0.4152 


21.85 


0.43 


29 


-19.1 


44 ± It 


0.3 


E1512-201268 11 


0.3412 


21.77 


0.39 


29 


-18.6 


qq + 45 
00 - 84 


1.4 


E1512-202096 11 


0.4252 


20.33 


0.66 


23 


-20.5 


Z ° - 33 


2.0 


E1512-201125 a 


0.3823 


21.69 


0.74 


15 


-18.9 


59 + 79 
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-0 
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+0 


.40 


-0 


20 


+0 


.40 


-0 


10 


+0 


.20 


-0 


10 


+0 


.20 


-0 


.10 


+0 


50 


-0 


.40 


+0 


03 


-0 


.00 


+0 


.20 


-0 


.10 


+0 


.40 


-0 


50 


+0 


•Ml 


-0 


.00 


+0 
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-0 


.10 


+0 


.10 


-0 


.10 


+0 


.10 


-0 


.10 


+0 


30 


-0 


.30 


+0 


.10 


-0 


.10 


+0 


30 


-0 


.20 


+0 


20 


-0 


.30 


+0 


.40 


-0 


.30 



a Kinematically anomalous galaxy 
n Kinematically normal galaxy 
c Serendipitously observed cluster galaxy 

Note. — (1) Galaxy identification code. First part is the CNOC cluster name and the second part is the CNOC 
PPP number; (2) CNOC galaxy redshift; (3) CNOC observed Gunn r magnitude; (4) CNOC observed Gunn g — r 
color; (5) Rest-frame [Oil] equivalent width. Obtained by dividing the observed width by (1+z). Error ~ 10%; (6) 
k-corrected absolute B-band magnitude (Ho = 75 km s _1 Mpc -1 , qo = 0.5). Error = ±0.1 mag; (7) Rest-frame 
projected [Oil] disk rotation velocity with 68% confidence interval; (8) [Oil] exponential disk scale length with 68% 
confidence interval (h = H o /(100 km s _1 Mpc -1 ), q = 0.5). Galaxies with no superscript were rejected because 
cosmic rays hit the [Oil] line directly, or simulations showed that the [Oil] flux was too low to derive reliable 
kinematical parameters. 
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Fig. 1. — Kinematical evidence for luminosity evolu- 
tion at intermediate redshifts. The locus of the local 
H Q rotation velocity— B band Tully-Fisher relation for 
all morphological types as defined by data taken from 
Mathewson et al. (1992) is used as a reference. Solid 
circles are the V rot sin i 's versus rest-frame B mag- 
nitude for the kinematically normal galaxies in our 
sample. The upper long dashed line is an unweighted 
linear fit to all the local morphological types. This lin- 
ear fit was then shifted by AMb = — 1.0 mag (middle 
dashed line) and AMb = —2.0 mag (lower dashed 
line) to represent various degrees of luminosity evolu- 
tion. 
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